Comparative study of strong coupling theories of a trapped Fermi gas at unitarity 
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We present a systematic comparison of the most recent thermodynamic measurements of a trapped 
Fermi gas at unitarity with predictions from strong coupling theories and quantum Monte Carlo 
(MC) simulations. The accuracy of the experimental data, of the order of a few percent, allows a 
precise test of different many-body approaches. We find that a Nozieres and Schmitt-Rink treatment 
of fluctuations is in excellent agreement with the experimental data and available MC calculations 
at unitarity. 
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The theory of strongly interacting fermions is of great 
interest. Interacting fermions are involved in some of the 
most important unanswered questions in condensed mat- 
ter physics, nuclear physics, astrophysics and cosmology. 
Though weakly-interacting fermions are well understood, 
new approaches are required to treat strong interactions. 
In these cases, one encounters a "strongly correlated" pic- 
ture which occurs in many fundamental systems ranging 
from strongly interacting electrons to quarks. 

The main theoretical difficulty lies in the absence of 
any small coupling parameter in the strongly interacting 
regime, which is crucial for estimating the errors of ap- 
proximate approaches. Although there are numerous ef- 
forts to develop strong-coupling perturbation theories of 
interacting fermions, notably the many-body T-matrix 
fluctuation theories [J, S, S, i, S, @, 0, S, i, their 
accuracy is not well-understood. Quantum Monte Carlo 
(QMC) simulations are also less helpful than one would 
like, due to the sign problem for fermions lj| or, in the 
case of lattice calculations l^, 3 , the need for extrapo- 
lation to the zero fllling factor limit. 

Recent developments in ultracold atomic Fermi gases 
near a Feshbach resonance with widely tunable inter- 
action strength, densities, and temperatures have pro- 
vided a unique opportunity to quantitatively test differ- 
ent strong-coupling theories 13, 15, 1^, 17, l3|- In these 
systems, when tuned to have an inflnite s-wave scattering 
length - the unitarity limit - a simple universal thermody- 
namic behavior emerges l^, 23|- Due to the pioneering 
efforts of many experimentalists, the accuracy of ther- 
modynamic measurements at unitarity has improved sig- 
niflcantly. A breakthrough occurred in early 2007, when 
both energy and entropy in trapped Fermi gases were 
measured without invoking any specific theoretical model 
[l^ . This milestone experiment, arguably the most ac- 
curate measurement in cold atoms, has an accuracy at 
the level of a few percent. 

In this Rapid Communication, using experimental data 
as a benchmark, we present an unbiased test of several 
strong coupHng theories that are commonly used in the 



literature, including QMC simulations. From this com- 
parison, we show that the simplest theory which incor- 
porates pairing fiuctuations appears to be quantitatively 
accurate at unitarity. This is the T-matrix approxima- 
tion pioneered by Nozieres and Schmitt-Rink (NSR) [l| 
and others 0, [H| , as recently appHed to trapped gases in 
the below threshold superfiuid regime We find it de- 
scribes the observed thermodynamics extremely well at 
all temperatures at unitarity, except in regions very near 
the superfiuid transition temperature T^. 

The comparisons show that the simple NSR approxi- 
mation gives excellent results. This appears to be related 
to the important symmetry property of scale-invariance 
2ll |. which is a necessary feature of any exact theory 



at unitarity, and is shared by the NSR approach. Our 
comparative results should therefore be useful in devel- 
oping new theoretical approaches for strong interacting 
fermions, and are relevant to many fields of physics. In 
particular, our results might shed light on the applicabil- 
ity of different T-matrix approximations to high- Tc su- 
perconductors and neutron stars, which are of interest to 
the condensed matter and astrophysics communities. 

The strong coupHng theories that we compare include 
several T-matrix fiuctuation and QMC theories. We re- 
fer to Refs. (l^ . [l3| for a detailed description of QMC 
methods, and briefiy review different T-matrix theories. 
These involve an infinite set of diagrams — the ladder 
sum in the particle-particle channel. It is generally ac- 
cepted that this ladder sum is necessary for taking into 
account strong pair fiuctuations in the strongly interact- 
ing regime, since it is the leading class of all sets of dia- 
grams j22| . 

The diagrammatic structure of different T-matrix the- 
ories may be clarified above Tc for a single-channel 
model [7], where one writes for the T-matrix, t{Q) = 
U/[l + UxiQ)]- Here and throughout, Q = (q, 
K = iKiujm), while = m/(47rft2a) - ^^m/ft^k^ 
is the bare contact interaction expressed in terms of the 
s-wave scattering length. We use J2k — ^BTj^m^k^ 
where q and k are wave vectors, while and tOm are 



2 



bosonic and fermionic Matsubara frequencies. 

Different T-matrix fluctuation theories differ in their 
choice of the particle-particle propagator (j^ . 

x{Q)=Y.j,G'-i^)GpiQ-K), (1) 

and the associated self-energy, S(i^) = J2q t{Q)G~f{Q ~ 
K). The subscripts a, /3, and 7 in the above equa- 
tions can either be set to "0", indicating a non-interacting 
Green's function Go{K) — — ft^k^/2m + fj], or 

be absent, indicating a fully dressed interacting Green's 
function. In these cases a Dyson equation, G{K) = 
Gn{K)/[l — Go{K)Ti{K)\, is required to self-consistently 
determine G and S. The only free parameter, the 
chemical potential ^, is fixed by the number equation, 
N^2Y.^G{K). 

By taking different combinations of a, (} and 7, there 
are six distinct choices of the T-matrix approximation, 
for which a nomenclature of {GaGp)G-y will be used. As 
noted earlier, there is no known a priori theoretical justi- 
fication for which is the most appropriate. While having 
the same diagrammatic structure, the T-matrix approx- 
imations we use above and below are computationally 
different, owing to the use of different Gq (or G). Below 
Tc, these Green's functions are 2x2 matrices. 

The simplest choice, {GaGo)Go, was pioneered by NSR 
above Tc using the thermodynamic potential ff)], with a 
truncated Dyson equation for G, i.e., G = Gq + GqSGo- 
This theory was extended to the broken-symmetry su- 
perfiuid phase by several authors 0, d, M, 0| ) using the 
mean-field 2x2 matrix BCS Green's function as "Go". In 
Ref. [1| , it was shown that the resulting ground state en- 
ergy is in excellent agreement with the zero-temperature 
QMC calculation for all interaction strengths. The NSR 
approximation is the simplest scheme that includes the 
effects of particle-particle fiuctuations. It does not at- 
tempt to be self-consistent. In the other extreme, one 
may consider a (GG)G approximation, with a fully self- 
consistent propagator. This was investigated in detail 
by Haussmann et al, both above and below Tc [lo| . 
Below Tc an ad-hoc renormalization of the interaction 
strength is required to obtain a gapless phonon spectrum. 

We also consider an intermediate scheme having an 
asymmetric form for the particle-particle propagator, 
i.e., {GGq)Gq. This has been discussed in a series of 
papers by Levin and co-workers based on the as- 
sumption that this treatment of fluctuations is consistent 
with the simpler BCS theory at low temperatures. Al- 
though the theory has been explored numerically to some 
extent ji^l, a complete numerical solution has not been 
implemented previously. A simplified version Q of the 
(GGo)Go fiuctuation theory was introduced based on a 
decomposition of the T-matrix t{Q) in terms of a con- 
densate part and a pseudogap part. In this Letter, we 
refer to this approach as the "pseudogap model" and will 
include it in our comparative study. 
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Figure 1: (color online) Chemical potential of a uniform Fermi 
gas at unitarity as a function of reduced temperature T/Tf, 
where Tf = eF/ks. The lines plotted are the results of NSR 
(solid lines), (GGo)Go (dashed line), {GG)G (dotted line), 
and pesudogap model (dot-dashed line). These predictions 
are compared with lattice QMC simulations (symbols). 

Other strong coupling theories with an artificial small 
parameter have been developed recently, including an e 
expansion around the critical dimension [H] and a 
expansion for a 2Ar-component gas ji^. These field- 
theoretic approaches provide very useful but so far only 
qualitative information about universal thermodynam- 
ics valid at unitarity. In the Boltzmann regime at high 
temperatures, not explored experimentally so far, it is 
possible to make a virial expansion in terms of fugacity 
[27|. We have verified that the three T-matrix schemes 
we study here do correctly include the dominant second- 
order virial contribution in the high temperature region. 

Fig. [T] compares the temperature dependence of the 
chemical potential at unitarity, calculated from different 
T-matrix schemes and lattice QMC simulations. The T- 
matrix approximations above Tc have been solved using 
an adaptive step Fourier transform method. Below Tc, 
the NSR and (GG)G results are from Refs. [1| and [lo| . 
respectively. The (GGo)Go approximation below Tc has 
not been worked out yet. The QMC results are taken 
from Refs. and [13|. However, these lattice calcula- 
tions may have systematic errors due to an extrapolation 
to the zero fiUing factor Hmit which is necessary to have 
a well-defined continuum theory. Nonetheless, the three 
T-matrix calculations agree well with the lattice QMC 
simulations. On the other hand, the prediction of the 
pesudogap model above Tc differs substantially from the 
(GGo)Go results, for which it is an approximation. The 
pseudogap model omits important features of the full 
(GGo)Go theory, due to the "condensate"+"pseudogap" 
decomposition of the T-matrix. 

The determination of energy and entropy is a sub- 
tle problem. It is known that universal thermodynam- 
ics at unitarity implies a rigorous scaling relation [l9| . 
P — —Vl = 2E/'i, which relates the pressure (or thermo- 
dynamic potential) and the energy density for a unitarity 
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Figure 2: (color online) Temperature dependence of the en- 
ergy (upper panel) and of the entropy (lower panel) of a uni- 
form Fermi gas at unitarity, obtained from different T-matrix 
approximations and QMC simulations as indicated. 



gas in the same way as for an ideal, non-interacting quan- 
tum gas, although the energy densities are quite different. 
Apart from the {GG)G scheme above-Tc and the NSR 
approach (in both regimes), strong coupling theories in 
general do not satisfy this essential scaling relation. The 
violation is typically at the level of a few percent, compa- 
rable to the accuracy of the experimental data we used. 
For quantitative purposes, we calculate the thermody- 
namic potential from the chemical potential, using 

n{ii,T = const) = - r n (/i') d^i' + n{^io, T) (2) 

J Mo 

at a given temperature. Here, the lower bound of the 
integral /io is sufficiently small so that ^{pLo^T) can be 
obtained accurately from a high temperature virial ex- 
pansion • The energy and entropy can then be cal- 
culated from the rigorous scaling relations, E = — 3r2/2, 
and S = (-5rj/2 - nN)/T, valid at unitarity. 

The energy and entropy obtained in this manner are 
given in Fig. [21 and compared to the predictions of QMC 
calculations. There is a reasonable agreement between 
T-matrix theories and the lattice QMC simulations. For 
the energy, we also show the path-integral Monte Carlo 



results of Akkineni et al. for the continuum model 
At temperatures above 0.2Tf, the energy lies systemat- 
ically below that of all the T-matrix theories. This is 
probably due to the use of a finite effective range tq for 



the interaction [ll|, i.e., kpro ~ 0.3. Compared to the 



QMC results, the pseudogap model appears to provide 
the least accurate predictions for energy and entropy. At 
low temperatures it predicts a T^/^ dependence for the 
entropy, which is characteristic of a non-interacting ideal 
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Figure 3: (color online) Theoretically predicted universal 
thermodynamics in comparison with experimental data [l8|]. 



Bose condensed gas. In contrast, the T-matrix entropies 
follow a scaling law, arising from the Bogoliubov- 
Anderson phonon modes flo| . 

We now compare theoretical predictions with exper- 
imental data [18|. A strongly interacting Fermi gas of 
N = 1.3(2) X 10^ lithium atoms is prepared in a Gaussian 
trap y(r) = Vb{l - exp[-m{uj^^p^ + lu'^^z'^)/{2Vo)]} with 
Vb ~ lOEp at a magnetic field B = 840 G, slightly above 
the resonance position Bq = 834 G. The large but finite 
interaction, kpa — —20.0, leads to an approximately 1% 
correction that is not accounted for experimentally. The 
energy is determined in a model-independent way from 



the mean square radius (z^) of the strongly interacting 
cloud, according to the virial theorem, which states that 
the potential energy (cx (z^)) of the gas is a half of its 
total energy. The entropy of the gas is calibrated again 
from the cloud size, but after an adiabatic sweep to a 
weakly interacting gas with kpa = —0.75, using a precise 
theory at weak coupling. We refer to Refs. and jlBl 
for further details. To determine the energy and entropy 
theoretically, we apply the local density approximation 
(LDA) by assuming that the system can be treated as 
locally uniform, with a position-dependent local chem- 
ical potential fihomin (r) ,T/TF{n)] = fJ. — V{r), where 
Tpin) is the local Fermi temperature. From this condi- 
tion, the density profile is obtained, and the total energy 
and entropy are calculated. 

Fig. [3] shows the interaction energy vs. entropy in a 
harmonic trap as predicted by the strong-coupling theo- 
ries in comparison with experimental data. All results of 
perturbation theories, except that of the NSR approach, 
were not reported previously to our knowledge. The en- 
ergy is expressed in units of the Fermi energy at zero 
temperature in a harmonic trap: Ep — (3A^w^a'z)^/'^ = 
ksTp. To emphasize the effects of interactions, we have 
subtracted the ideal gas result Ejg- No adjustable pa- 
rameters have been used theoretically or experimentally. 
This comparison is therefore an unbiased test of how well 
T-matrix theories agree with experiment [l8| . 
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The difference between different T-matrix schemes, 
mostly of the order 0.05NEp, is relatively small and thus 
nearly indistinguishable in the plot of total energy and 
entropy. Despite this, the extraordinary precision of the 
measurements is able to discriminate between these theo- 
ries in the interaction energy, as given in Fig. [3l The NSR 
approach is seen to give the best fit to the experimen- 
tal data below Tc (corresponding to Sc « 2.3NkB) and 
above T = O.STp (corresponding to S > S-dNks)- This 
indicates that the simplest non-self-consistent T-matrix 
approximation captures the essential physics of strong 
pair fiuctuations at both low (superfiuid) and high (nor- 
mal) temperatures. Around Tc, the experimental data 
shows evidence of what could be a first-order superfiuid 
transition. However, due to "critical slowing-down", sys- 
tematic experimental errors cannot be ruled out in this 
regime, if the magnetic field sweep is not quite adiabatic. 

In the temperature region just above Tc, the NSR ap- 
proach presumably does not fully capture the full effect of 
fiuctuations, compared to the self-consistent {GG)G the- 
ory above Tc. At the transition, from the experimental 
data one may determine experimentally a critical entropy 
Sc/N ~ 2.3kB and a critical energy Ec/N ~ 0.86£'f in 
a trap. The critical temperature Tc/Tp in the case of 
a trap is difficult to determine, due to the unknown re- 
lation S{T). The theoretical predictions are 0.29 [NSR], 
0.21 [(GG'o)Go and (GG)G], and 0.27 [pseudogap model]. 

In a further comparison, we include in Fig. [3] a recent 
QMC calculation (thin solid line) of trapped Fermi gases 
[231 ■ There is a noticeable systematic difference between 
the QMC and experimental data at high entropy, but this 
is due to the improper use of an ideal gas approximation 
in the QMC estimates for large T. It is clear that the uni- 
tarity gas remains strongly interacting even close to the 
degenerate temperature {i.e., S ~ ^Nks)- Thus, a virial 
expansion of the equation of state up to the second order 
should be applied. Among all pair fiuctuating theories. 
Fig. 3 shows that the pseudogap approximation gives 
poor agreement with thermodynamic data, though it is 
better than BCS mean-field theory - which completely 
ignores the pairing fiuctuations. Therefore, the pseudo- 
gap model does not describe the strong fiuctuations at 
unitarity as well as the full (GGo)Go theory. 

In conclusion, the accurate thermodynamic measure- 
ments at Duke University have allowed us to perform 
a test of strong-coupling T-matrix theories at unitarity. 
The simplest NSR approximation for the particle-particle 
T-matrix is found to give the best quantitative descrip- 
tion. Further work is needed to understand the reason 
for this, but we conjecture that it is related to scale- 
invariance symmetry in the unitarity limit. We conclude 
that near the BCS-BEC crossover region, the NSR ap- 
proximation is surprisingly useful. 
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